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Abstract 


Vertical  radiation  patterns  of  vertical  quarter -wave  monopoles, 
and  of  vertical  and  horizontal  half-wave  dipoles  located  at  various 
heights  up  to  twenty  wavelengths  over  several  typical  types  of  ground 
are  shown.  Vertical  antennas  have  the  best  low-angle  patterns  at  low 
heights  and  at  medium  heights  over  land.  F or  low  and  medium  antenna 
heights,  only  vertically  polarized  antennas  over  water  radiate  appre¬ 
ciable  energy  at  zero  or  negative  elevation  angles. 

The  propagation  factor  used  to  calculate  the  antenna  patterns  can 
be  accurately  evaluated  with  the  simple  two-ray  theory  for  elevation 
angles  greater  than  one  or  two  degrees.  For  smaller,  and  negative, 
elevation  angles  the  full  wave  theory  for  diffraction  by  a  sphere  must 
be  used. 


Keywords:  Antenna  patterns,  HF  Antennas,  low-angle  radiation, 
propagation  factor. 
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Linear  High  Frequency  Antennas  Over  a 
Finitely  Conducting  Spherical  Earth 

Leslie  A.  Berry  and  Mary  E.  Chrisman 
1.  Introduction 

Recently  there  has  been  increased  interest  in  using  low 
elevation  angle  paths  for  HF  ionospheric  communications.  Utlaut 
[1961]  measured  HF  signals  propagated  over  6800  km  and  8400  km 
paths.  On  the  average,  more  energy  arrived  at  low  angles  (<5°)  than 
at  high  angles  (~150),  even  though  the  transmitting  antenna  favored  high 
angles.  F or  a  given  signal  threshold,  Utlaut  found  that  low-angle 
reception  provided  service  a  greater  percentage  of  the  time.  Other 
reports  by  various  workers  suggesting  the  same  result  are  listed  and 
discussed  by  Utlaut. 

Because  of  the  possible  usefulness  of  low-angle  propagation, 
the  calculation  of  low-angle  radiation  patterns  of  HF  antennas  is 
desirable.  Commonly,  antenna  patterns  are  computed  assuming  that 
the  earth  is  flat  and  smooth  (and,  often,  perfectly  conducting).  Both 
of  these  assumptions  are  questionable  at  low  elevation  angles,  and  are 
certainly  invalid  at  zero  or  negative  elevation  angles.  In  this  report, 
the  model  of  a  flat,  smooth  earth  is  replaced  by  a  spherical,  smooth 
earth,  and  some  consideration  is  given  to  the  effect  of  a  large  surface 
irregularity. 

The  necessary  solution  of  Maxwell1  s  equations  has  existed 
for  some  time,  since  the  basic  problem  is  just  that  of  diffraction 
by  a  conducting  sphere.  The  first  useful  solution  for  a  large  sphere 


was  given  by  Watson  [1918].  Bremmer  [1949]*  and  Fock  [1965]^  made 
significant  progress  in  providing  useful  approximations  for  numerical 
work.  The  functions  now  most  frequently  used  to  express  the  solution 
were  defined  by  Fock  although  equivalent  solutions  were  given  by  Pryce 
[l953]in  a  different  notation.  Norton  [  1 941  ]  published  graphs  which 
made  simple  engineering  calculations  of  the  ground  wave  possible.  Wait 
and  Conda  [1958]  gave  numerical  values  of  the  current  induced  on  a  curved 
conductor  by  an  incident  plane  wave  --  corresponding  to  a  short  dipole  on 
the  surface  receiving  a  downcoming  sky  wave.  In  1959,  they  considered 
the  far  field  diffraction  by  a  cylinder  of  an  incident  plane  wave  [Wait  and 
Conda,  1959].  Logan  [1959]  gives  an  excellent  historical  review  of 
previous  work,  lists  over  80  original  contributions  to  the  problem,  and 
shows  the  relations  between  the  functional  notations  used  by  the  various 
author  s . 

In  this  paper,  the  antenna  is  within  a  few  wavelengths  of  the  surface 
and  a  target  source  is  located  in  the  far  field.  The  vertical  radiation 
patterns  of  practical  antennas  are  found  by  integration  over  the  field  for 
a  current  element.  Example  calculations  are  made  for  vertical  and 
horizontal  half-wave  dipoles,  and  vertical  quarter -wave  monopol'js  for 
a  frequency  of  10  MHz. 

The  conclusions  are:  In  the  HF  band,  the  curvature  of  the  earth 
is  important  only  for  elevation  angles  less  than  one  or  two  degrees. 

Except  for  vertical  antennas  over  water  and  very  high  antennas,  little 
energy  is  received  at  these  angles,  so  the  curvature  of  the  earth  is  at 
best  a  second  order  effect.  As  expected,  the  best  low-angle  radiation  by 

1.  The  book  by  Bremmer  summarizes  and  extends  the  work  done  by 
Bremmer  and  Van  Der  Pol. 

2.  Fock's  research  was  done  in  the  1940s.  The  book  referenced  is 
an  excellent  English  translation  of  his  earlier  papers  on  diffraction 
theory. 


2 


low  antennas  is  by  vertical  polarization  over  sea  water.  The  low- 
angle  patterns  of  low  horizontal  antennas  improve  as  the  ground 
conductivity  decreases.  For  high  antennas,  the  low-angle  patterns  of 
horizontal  antennas  are  as  good  as  those  of  vertical  antennas. 

2.  The  Propagation  Factor 

The  formulas  for  the  antenna  patterns  are  conveniently  written 
in  terms  of  the  ''propagation  factor",  which  is  the  ratio  of  the  electric 
field  at  a  point  to  the  field  which  the  same  source  would  excite  at  that 
point  if  all  boundaries  were  removed;  i.  e.  ,  the  ratio  of  the  actual  field 
to  the  free  space  field.  This  factor  has  also  been  called  the  reflection 
factor,  Laitinen  [1949],  the  cutback  factor  [Wait  and  Conda,  1958], 
and  the  attenuation  function  [Norton,  1936;  Furutsu,  1965]. 

2.  1  Flat  Earth,  Plane  Wave  Example 


As  a  simple  example,  consider  a  probe  antenna  located  at 

height  h  above  a  plane  earth  as  shown  in  figure  1.  A  plane  wave  E0 

is  incident  from  the  right  at  an  angle  A  with  the  horizontal.  When  the 

lower  ray  strikes  the  ground  some  of  its  energy  is  absorbed,  and  some 

is  reflected  with  a  phase  shift.  The  ratio  of  the  reflected  wave  to  the 

incident  wave  is  the  reflection  coefficient,  R,  i.e.  ,  E  _  =  RE0. 

reflected 

In  general,  R  is  a  complex  number.  Then  the  total  field  received  at 
A  is 


E  =  E0  +  e'lCPRE0  =  E0(l  +  e"lCPR), 


(1) 


where  cp  is  the  difference  in  phase  path  between  the  direc  t  and  reflected 
rays.  From  geometry, 

Cp  =  k  2h  sin  A,  (2) 
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where  k  =  2 rr/\  and  X  is  the  free  space  wavelength.  For  this  case 
-icp 

F(A,h)  =  1  +  e  R  is  the  propagation  factor. 

The  reflection  coefficient  R  is  a  function  of  the  angle  A,  the 
polarization  and  frequency  of  the  incident  wave,  and  the  electrical 
properties  of  the  reflecting  surface.  For  a  surface  with  relative 
dielectric  constant  e,  and  conductivity,  a  mho/m,  the  reflection 
coefficient  for  "vertical  polarization"  (electric  vector  in  the  plane 
of  incidence)  is  [Norton,  1936] 


where 


k2 


18J(103) 

fMHz 
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and  |i0  —  4tt(10  )  is  the  permeability  of  free  space,  c  is  the  speed  of  light, 

oj  =  2nf  is  the  angular  frequency,  and  is  the  frequency  in  megahertz. 

For  horizontal  polarisation  (electric  vector  perpendicular  to  the  plane  of 


Some  characteristics  of  ground  reflection  coefficients  will  be  dis¬ 
cussed  later.  The  most  important  thing  for  HF  is  that  if  |ke  |  »  k 
(high  conductivity),  R^  —  -  1 ,  and  except  for  very  small  A,  R^  ^  1. 
If  A  =  0.  1^  =  1^=  -  1. 
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2.  2  Spherical  Earth,  Vertical  Component  in  the 
Plane  of  Incidence 

We  will  now  develop  formulas  for  the  propagation  factor  for  a 
sky  wave  incident  on  a  spherical  earth.  The  geometry  of  the  spherical 
model  is  shown  in  figure  2.  The  receiving  antenna  is  an  elementary 
dipole  at  A,  hA  meters  above  a  smooth,  homogeneous,  spherical 
earth  of  radius  a  whose  center  is  the  origin  of  the  spherical  coordinate 
system  (r,  0,  cp).  To  model  waves  reflected  from  the  ionosphere,  a 
target  source  T  is  located  at  ionospheric  reflection  height  h,^,  a 
distance  D  from  A.  Throughout  this  report  we  assume  h  >  h  .  The 

X  A 

direct  ray  TA  makes  an  angle  A  with  the  local  horizontal  at  A  --  i.  e.  , 

A  is  the  elevation  angle.  The  reflected  ray  travels  the  path.  Di  +  D2 . 

The  angular  distance  between  A  and  T  is  0  =  —  . 

8l 

For  large  elevation  angles,  the  flat  earth  approximation  developed 
above  may  be  used  by  defining  the  equivalent  height  fi  =  h  +  D  sin  A. 

X  A 

Let  the  target  source  at  T  be  a  vertically  polarized  elementary  dipole 
with  current  moment  I  dh,  where  dh  is  along  the  radius  of  the  sphere. 
Then  the  radial  electric  field  at  A  (assuming  time  dependence 
exp(iuut)  which  is  not  written)  is  [Feldman,  1933;  Laitinen,  1949] 


E  ^  -  30ik  ---—I  dh  cossa(i  +  R  e-lk<Di  +  Qa  ~Dh  <  (5) 

r  D  V  V  y 

This  equation  is  (-cos  A)  times  the  one  usually  given  because  it  gives 
the  E^  component  for  a  coordinate  system  originating  at  the  earth*  s 

center  instead  of  the  E  component  for  a  coordinate  system  centered 

0 

on  the  antenna.  For  large  D  and  h^  »  h^  f 

Dx+Pa  -D  as  ZhA  sin  A  f  (6) 

which  reduces  (5)  to  (1). 


> 
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As  A  approaches  zero,  the  notion  of  distinct  direct  and  reflected 

rays  breaks  down.  Waves  coming  from  T  graze  the  horizon  and  travel 

some  distance  near  the  curved  ground  to  get  to  A.  In  fact  there  is  a 

continuum  of  waves  which  reach  A  in  this  manner,  and  they  must  be 

summed  to  obtain  the  resultant  field.  Thus  w'e  must  integrate  over 

this  continuum  and  use  wave  functions  appropriate  to  a  spherical 

geometry.  Exactly  how  this  should  be  done  would  not  be  clear  without 

referring  to  the  rigorous  direct  solutions  given  by,  for  example,  Fock 

[1965]  and  Wait  [1962].  Since  the  intent  here  is  to  stress  the  physical 

interpretation,  the  appropriate  wave  functions  and  form  of  the  integral 

are  taken  from  these  references  without  derivation. 

Wave  functions  appropriate  to  the  problem  are  the  Airy  functions 

W  (t)  used  by  Fock  [1965],  Wait  [1962],  and  others.  The  functions  are 
n 

linearly  independent  solutions  of  the  differential  equation, 

W'  (t)  =  tw  (t)  .  (7) 

n  n  • 

for  n  =  1  or  2,  and  for  -t  »  1, 

WiW^t-t)  4  exp  [T  i(z  +  tt/4)^  ^ 

and  \  (8) 

Wi(t)  —  (-t)4e:  p  [^i(z  -  n/4)].  I 

3  j 

where  z  =  §•  (-t)  2,  Equation  (8)  shows  that  the  functions  act  like 

plane  waves  for  large  (-t)  and  W1  (t)  is  an  upgoing  wave  and  W2  (t)  a 
downgoing  wave. 

The  relation  of  the  variable  t  to  the  geometry  is 

t  ss  -v2  sin2  A  ,  A  >  0  ,  (9) 
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1/ 

where  v  =  (ka/2)  3  .  (10) 

The  so-called  natural  units  of  distance  and  height  are  [Pryce, 

1953] 


x  =  v9 ,  and  y  =  kh/v.  (11) 

Then  the  integral  which  sums  the  continuum  of  waves  reaching  A 
from  T  at  small  angles  A  is  (compare  Wait  and  Conda  [1958]  and 
Fock  [1965] ) 

,  ikd 

E  2=e  '  KU- -  I  ,  and  (12a) 

r  « /  sin  0  r 


Vs 


Wi(BA)  W2(t) 


r*  wt/  f  n  a  2'  '  A  1 

Vf/  l1  +  27-)  wTioS)11'-  (12b| 

1  A 


In  (12),  K  =  12  vVk/a3"  I  dh  , 


u  =(4)U/^i 


(13) 

(14) 


where  r  =  a  +  h  ,  b  =  a  +  h  ,  and  T  is  a  contour  enclosing  the  poles 

A  J. 

of  the  integrand  in  the  lower  half  plane.  For  discussion  here,  T  can 
be  taken  along  the  real  axis  from  ®  to 

.  a  ./3 

Bn=(irr)  (‘  •  Yn)>  n  =  A  or  T  , 
n 

—  (t  -  y  ),  since  h  «a.  Most  authors  use  the 
n  n 

latter  form,  and  it  will  be  used  for  the  rest  of  this  section,  but  the 
more  accurate  form  was  necessary  for  some  of  the  numerical 
calculations. 
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Finally,  R^  is  the  spherical  reflection  coefficient, 

A  W2'(t)/w2(t)  -  qv 

Rv  =  "  wTTtjTwHtj  -  qv  ’ 


where 


k 

q,r  =-ivr 


Equation  (12)  has  only  a  superficial  resemblance  to  (5),  so  the 
first  thing  to  do  is  to  show  that  it  indeed  reduces  to  (5)  for  large  A#  This 
demonstration  will  also  show  the  roles  of  the  various  factors  in  (12). 
Referring  to  (9),  if  A  is  large,  (-t)  »1,  since  v  >  10  at  high  frequencies. 
So  substitute  (8)  into  (12)  and  simplify: 


{i+4~s)5/s 


-f  — - ,  +  S>Rv)dt. 

r  -«  (YA-t)  ^(Y  t)  V 


where  Q(t)  =  xt  +  §-  (  Y_  -  t)  ^  -  \  (  Y  -  t )  ^  . 


This  integral  can  be  evaluated  by  the  saddle  point  method  [Morse  and 
Feshbach,  1953].  The  result  is 


where 


S>  =  f-  ( Ya-  to)3/s  -t-  (-to)  /s,  and 
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to  is  the  solution  of  0/(to)  =  0,  or 


2x2(yA+yT)  ■  (yA'yT^  '  x4 


Then 


E  .Z^iKT  J TTe-i(kdW(1o))A+e-iq)&  > 

r  rrr~7r  x  V  yy 


sin  9 


which  now  has  a  form  much  like  (5). 

First,  we  want  to  verify  that  (20)  is  equivalent  to  (9)  to  justify 
the  transformation  of  variables.  Substitute  (11)  into  (20)  and  simplify 
to  get 


.  r<VV  ,  t  hA; 

tb  =  -v  {— JS - - J 


s2 

j  ,  if  d2  «  4a(hT-hA)  , 


which  corresponds  to  A  being  somewhat  greater  than  0. 

From  the  geometry  of  figure  2,  if  h  «  h  (the  antenna  is 

A  i 

much  lower  than  the  ionospheric  reflection  height), 


sin  A  =  1  -  cos2  A  —  1  -  tts  — 


d2  D2  -d2 


D  D2 


Since  D2  2=  d2  +  (hT-hA)2  » 


•  2  A 
sm  A 


(hT"hA)  (  r  hT'hA  ^  (hT'tlA) 

— —  -  v d — )  J  ~  “? - 


if  A  is  not  too  large. 
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Substitution  of  (23)  into  (22)  yields  (9). 
Now,  using  (8)  in  (15)  9 


i 


or,  using  ( 1 6)  and  (9), 


A 


Ir 

/  f  k  'f 

sinA-^V 

'  1  -  Vk J 

sin  A  +  —V 
ke 

(24) 


A 

Thus  Ry  —  Ry  if  cos  A  —  1,  which  is  true  for  small  A.  This  justifies 

A 

identifying  R^.  as  the  reflection  coefficient  and  shows  that  a  more 
accurate  definition  of  q^.  would  include  the  cos2  A  as  a  factor  of 
But  since  the  integral  is  to  be  used  only  for  very  small  A  and  since 
(k/kg  )2  is  usually  quite  small,  it  does  not  seem  necessary. 

Next  consider  the  phase  terra  kd  +  Q(to)  in  (21).  By  putting  (11) 
and  (20)  into  (18)  and  grinding  through  some  tedious  algebra,  we  get 


n{to)  =  k( 


d(hA+hT) 

2a 


+ 


<W 

2d 


(25) 


By  the  law  of  cosines,  =  r2  +  b2  -  2rb  cos  0.  Using  Taylor*  s 
series  for  cos  9,  after  working  out  the  algebra, 


d(h  A+h  )  (h  -li  )2  j3 


D  J-  A  T  |  A  T 


2a 


2d 


24a2 


+  terms 


<D‘- 


Thus,  to  a  close  approximation,  0(to)  =  k(D-d)  and 
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(26) 


-i(kd-Q(tb))  -ikD 
e  =  e 


Similarly,  cp  -  cp. 

Finally,  consider  the  coefficient: 


-2iKU 
*]  sin  0 


-24i  v®7k/a3  Jn/xQ  -  Sl°  ^ 

Jd/a. 


I  dh 


■  ■ 3oi  1 0  ■  i)s  /a  1 611  (z7) 


d  A2 

for  6  =  —  «1  and  A  «  1.  But  for  small  A,  d  —  D  cos  A  and  1  —  cos  A, 

a  2 

so  the  coefficient  is  approximately 


cos2 A  ^ 

-3  Oik  — —  I  dh  , 


This  reduces  (21)  and  hence  (12)  to  (5)  for  the  range  of  A  in 
which  both  are  valid.  It  may  seem  contradictory  that  we  have  required 
the  elevation  angle  A  to  be  both  "not  too  small"  and  "not  too  large". 
This  is  because  we  are  trying  to  show  that  two  formulas,  one  for 
large  A  and  one  for  small  A,  have  an  overlap  region  in  which  both 
give  the  same  answer.  Careful  examination  of  the  inequalities  above 
show  that  they  require  A  to  lie  between  about  2°  and  10°.  Numerical 
calculation  of  (5)  and  (12)  has  shown  that  they  do  indeed  give  the  same 
result  in  this  region. 

Comparison  of  (19)  and  (26)  shows  that  if  I  is  calculated 
directly  then  the  propagation  factor  is  given  by 
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Fr(A,hA)  =  -  £  JZ  e 
TTl 


-ik(D-d) 


(28) 


Reviewing  the  results  of  this  section  so  far,  it  has  been  shown 
that  the  simple  two  ray  formula  for  the  propagation  factor  can  be 
derived  from  the  rigorous  full  wave  solution  for  sufficiently  large  ele¬ 
vation  angles  in  accordance  with  physical  intuition.  For  very  small  and 
negative  elevation  angles,  a  continuum  of  waves  must  be  summed  with 
the  integral  shown  in  (12).  The  integral  can  be  evaluated  numerically 
on  a  computer  by  swinging  the  ends  of  the  contour  down  into  the  lower 
half  plane  to  get  convergence,  being  careful  not  to  cross  any  poles  of 
the  integrand. 

The  integral  in  (12)  can  also  be  evaluated  with  the  calculus  of 
residues  for  negative  take-off  angles.  Rewrite  (12)  as 


(l+7”g)  -ixt 

* 2i  I  w;it)7wV(t)  -;v  h,<,ti  «■<  V dt 


(29) 


where  Hi(y_)  = 


W1(t-yT) 


T'  Wi(t) 


,  and 


H3(yA)  = 


Wa(t-yA)  [w;<t)  -  qvVM t)]  -  W1(t-yA)  [W2'(t)  -  qyW2(t)] 


This  form  has  two  advantages:  first,  the  H^y)  are  height  gain  functions 

such  that  if  h  =  0,  y  =  0  and  H  (0)  =  1.  This  is  obvious  for  Hx  (y). 

n 

For  H2(y)  it  is  shown  by  using  the  Wronskian  [Wait,  1962], 


W/WW^t)  -  Wa  (t)W2'(t)  =  2i.  (30) 

The  second  advantage  is  that  this  form  shows  that  the  integrand  has 
poles  at  the  points  tg  which  are  solutions  of 
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(31) 


Wl'(ts)  =  qVWl(ts)- 


Then  the  field  is  given  by 


Z  5  / 

CO  c  a/_  — ivt 

„  .  A  -3M/4  r  <‘*2?  '  *  ) 

E  a  -  4rre  KT -  ) 


^^*  =  0  ls  ’  q% 


Wi(t  )  Wi(t  ) 
s  s 


(32) 


which  is  the  well-known  series  for  the  ground  wave  (shown  in  this 
notation  by  Wait  [1962]  andFock  [1965]).  Equations  (30)  and  (31) 
were  used  to  simplify  Hg .  In  this  case,  the  propagation  factor  is 


FjA.hJ  =  2rr  /£•  e"lk(D"d^  (residues). 

s 


fTTi 


(33) 


2.  3  Horizontal  Component  in  the  Plane  of  Incidence 

The  propagation  factor  for  the  horizontal  component  in  the 
plane  of  incidence  of  an  incoming  wave  is  [Laitinen,  1949] 


F0(A,h)  =  i  -  e‘lCPRv. 


(34) 


In  the  spherical  geometry  of  figure  2,  the  Eq  component  is  such  a 

wave.  The  integral  corresponding  to  (29)  for  EQ  is  (for  h_  >  h  ) 

u  T  a 

3/- 

!/3„e‘ikd  r  0+zW 


-ixt 


E.«  ■.“'“zngY V  _ 

v/sin  0  W/ (t)/Wx  (t)-q 


j. 


Hl(yT)H3(yA)  dt 


(35) 
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where 


Hs(yA)  = 


w i (t-yA)  [w; (t)-qyw,  (t)]  -  Wi' (t-yA)  [w2' (t)-qyW2  (t)] 


(36) 


By  evaluating  (35)  with  the  saddle  point  approximation  it  can  be  shown 
that  it  reduces  to  the  plane  wave  formula  just  as  in  the  previous  case. 
The  propagation  factor  is 


Fe<A’hA> 


vv^ 10 


(37) 


where  Iq  is  the  integral  of  (35).  The  denominator  of  (37)  is  zero  when 


d  =  ,/2a(hT-hA) 


(38) 


which  occurs  when  A  =  0.  But  there  is  no  infinity  in  the  field 

Eq  =  E  (free  space)  *  Fq  because  the  9 -component  of  E  (free  space) 

is  proportional  to  sin  A. 

2.  4  Horizontal  Component  Perpendicular  to 
Plane  of  Incidence 

The  term  "horizontal  polarization"  usually  indicates  that  the 
electric  vector  is  perpendicular  to  the  plane  of  incidence.  In  the 
coordinate  system  of  figure  2,  this  corresponds  to  the  E^  component. 
For  large  At  the  propagation  factor  is 

F^A.h)  =  (1  +  e"1C?Rh)  (39) 


where  is  given  by  (4).  Since  (39)  has  the  same  form  as  (1)  it  seems 
reasonable  that  the  analysis  of  section  2.  2  should  apply  here  just  by 


replacing  R^.  by  R^.  This  is  indeed  the  case,  and  the  integral  for 
has  the  same  form  as  (29)  except  that  q^.,  (16),  is  replaced  by 


\--lvryi-(£) 


(40) 


2.  5  Relation  to  Other  Problems  and  Special  Cases 
The  integral  in  (12)  and  (29)  and  the  residue  series  (32)  are  the 
same  as  those  given  by  Fock  11965  land  by  Wait  and  Conda  [1958]  except 
for  the  factor  ^1  +  ^  ”  1  for  small  t  (small  A).  Wait  [1962] 

showed  that  it  is  equivalent  to  the  ground  wave  series  of  Bremmer 
[1949].  If  the  receiving  antenna  is  put  on  the  surface (  h^  =  0/,  and  the 
source  is  raised  very  high,  the  formulas  correspond  to  those  given  for 
the  cut-back  factor  by  Wait  and  Conda  [1958].  If  both  antennas  are 
raised  very  high,  (12)  is  equivalent  to  the  diffraction  integral  denoted 

A 

G(X)  by  Wait  and  Conda  [1959]  and  also  given  by  Fock  [1965],  Thus 
the  solution  used  here  is  the  general  solution  for  electromagnetic  wave 
propagation  near  and  around  a  large,  finitely  conducting  sphere.  It 
contains  implicitly  both  the  space  wave  and  the  surface  wave  [Norton, 
1936],  although  it  has  not  been  written  to  display  them  separately. 

As  a  final  checK,  the  wave  tilt  (j=  —  j  [Norton,  1937]  will  be 
calculated.  Let  h  =  0,  so  that  the  fields  are  found  on  the  surface, 

A 

and  H3  =  Hi  =  1.  For  A  <  0,  only  a  diffracted  wave  is  present.  Then 
since  ( t  +  J  a:  1 ,  dividing  (35)  by  (29), 


1CL 


r 

J... 


dt 


1CL 


i  •  •  •  dt 


v 


k_  ! 

k2  J  1  k  k2  J 


(41) 


because  the  integrands  are  now  the* same.  The  result  on  the  right  is 
the  correct  value  tor  the  wave  tilt  [Norton,  1937], 
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2.  6  Calculation  and  Discussion  of  the 
Propagation  Factor 

The  propagation  factor  was  calculated  for  frequencies  from  1  to 
25  MHz,  for  the  types  of  ground  shown  in  table  1,  and  for  heights  h  = 

A 

.  05,  .1,  .25,  .375,  .  5,  .  75,  1,  2,  5,  and  10 


Table  1.  Ground  Constants  Used  in  this  Report 

Relative  dielectric 


Ground  type 

Conductivity,  mho/m 

constant 

sea  water 

5.0 

80 

good  ground 

r-H 

o 

• 

o 

10 

poor  ground  and 

sea  ice 

0.  001 

4 

polar  ice  cap 

0.  0001 

1 

fresh  water 

0.  002 

80 

The  values  for  sea  water,  fresh  wa'ter,  good  and  poor  ground 
are  fairly  standard;  that  for  sea  ice  is  an  average  taken  from  Wentworth 
and  Cohen  [1964],  The  constants  for  "polar  ice  cap"  are  speculative  -- 
they  at  least  represent  the  extreme  of  poor  conductivity. 


Printed  tables  and  magnetic  tapes  containing  the  data  are 
available.  To  illustrate  the  properties  of  the  propagation  factor,  the 
amplitude  has  been  plotted  for  f  =  1,  5,  15,  and  25  MHz  for  sea  water 
and  poor  ground  for  both  vertical  and  horizontal  polarization.  They 
will  be  discussed  later  in  this  section,  after  explanation  of  some 
details  of  the  calculations. 

The  integral  formulas  for  the  propagation  factor  for  small  ele¬ 
vation  angles  (29,  35)  apparently  depend  on  the  height  h^,  of  the  target 
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source,  an  undesirable  feature.  In  fact,  the  propagation  factor  depends 

on  li  only  when  T  is  below  the  horizon.  This  occurs  when  -A  >  2h.  /a. 

T  7  A 

Even  then,  the  dependence  on  h,^  is  very  slight  over  a  reasonable  range 
of  hT>  Still,  hT  must  be  given  some  value  in  the  formulas,  and  in 
this  report,  h,^  =  200  km. 

F or  a  number  of  test  cases  calculated  for  h,^  =  300  km  and  h^  = 

100  km,  the  difference  in  |f|  was  less  than  0.  5  dB  if  -1. 5°  ^  A  <  0.  5°, 
and  there  was  no  difference  if  A  >  0.  5°.  This  error  is  hardly  significant 
since  the  antenna  patterns  shown  in  section  3  are  all  more  than  20  dB 
below  an  isotropic  radiator  if  A  <  -1.  5°. 

The  decrease  in  the  refractive  index  of  air  with  height  causes 
some  bending  of  a  radio  wave.  The  amount  of  bending  varies  with  time, 
location,  and  elevation  angle*,  but  on  the  average  it  is  accounted  for  by 
using  an  equivalent  earth' s  radius  4/3  as  large  as  its  actual  mean  radius 
[Burrows,  1935].  This  standard  pxocedure  was  used  for  the  calculations 
of  the  propagation  factor. 

The  condition  for  use  of  the  saddle  point  approximation,  and 
hence  the  two-ray  model,  is  -to  —  (v  sin  A)2  »  1.  In  practice,  if 
-t0  >  4  (i.e.,  A  >  2l  at  1  MHz,  A>  1°  at  10  MHz),  the  two- ray,  plane 
wave  calculation  is  correct  to  at  least  3  significant  digits,  so  at  this 
vi  lue  the  computer  program  switches  to  numerical  integration  of  the 
integral.  When  A  <  0,  the  residue  series  (32)  is  used. 

Actually,  the  ranges  of  validity  of  the  three  methods  overlap, 
providing  a  check  on  the  logic  and  precision  of  each  branch  of  the 
program.  The  program  was  also  validated  by  comparing  with  calculations 
of  the  cutback  factor  by  Wait  and  Conda  [  1958]  and  Hyovalti  L  1  96  5 J ,  and 
with  the  ground  wave  calculations  of  Johler,  Walters,  and  Lilley  A  1  96  0  J , 
since  these  are  special  cases  of  the  general  solution  used  here.  The 
only  unresolved  difficulty  is  the  inability  to  make  a  smooth  transition 
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from  the  two-ray  solution  to  the  integral  for  h  =  20  \  at  1  MHz  (h  = 

A  A 

6000  m).  If  it  is  necessary  to  make  calculations  for  such  a  high  antenna, 

a  smooth  transition  could  be  interpolated. 

Figure  3.  1  shows  the  propagation  factor  for  the  vertical  and 

3 

horizontal  components  for  h  =  0.  05A.  over  sea  water  .  The  lines  are 

A 

labeled  with  the  frequency  in  MHz.  At  first  glance,  the  most  striking 

feature  is  that  there  is  almost  no  frequency  dependence  for  horizontal 

polarization,  and  very  little  for  vertical  polarization  above  A  =  10°. 

There  are  two  reasons  for  the  very  weak  frequency  dependence:  First, 

in  the  graphs  the  height  is  measured  in  wavelengths.  If  the  curves  were 

plotted  for  some  fixed  height  in  meters,  there  would  be  much  more 

dependence  on  frequency.  Second,  there  is  no  frequency  dependence 

for  a  perfect  conductor,  and  sea  water  is  a  good  conductor. 

Another  noticeable  feature  is  that  the  amplitude  of  the  propagation 

factor  is  much  greater  for  vertical  polarization  than  for  horizontal 

polarization.  For  good  conductivity,  Ry  —  1>  except  for  very  small 

elevation  angles,  and  R  ~  -1.  For  horizontal  polarization  the  reflected 

h 

ray  nearly  cancels  the  direct  ray,  while  for  vertical  polarization  the 
two  reinforce  each  other  and  1 1  +  Ry|  —  2. 

Finally,  note  that  the  vertically  polarized  wave  is  not  cut  off  at 
A  =  0  as  it  would  be  for  a  flat  model  of  the  earth,  but  is  merely  cut  back. 
For  negavive  elevation  angles,  the  amplitude  of  the  propagation  factor 
decreases  exponentially  with  frequency,  just  as  the  ground  wave  does. 

Figures  3.  2  and  3.  3,  for  heights  of  0.  1  A.  and  0.  25A.  are 
similar  to  figure  3.  1. 

The  propagation  factor  for  a  perfect  plane  reflector  is  [Laitinen, 

1949] 

3.  In  the  figure  legend,  HT  is  the  height  h  . 

A 
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**i  k  h  sin  A 

Fy  =  2  e  cos(khsinA),  and 


.  -lkhsmA  .  . 

F,  =  2  1  e  sm  (khsmA). 

n 


For  small  A,  the  first  null  of  F  is  at  sin  A  =  —  ,  and  the  first  null  of 
X.  ^  4x1 

F,  is  at  sin  A  =  —  .  If  h  =  X/2,  as  in  figure  3.  4,  F,_  should  have  a  null 
h  Zh  V 

near  A  =  30°  --  and  it  does.  For  the  same  height,  F  should  have  a  null 

h 

near  A  =  90°.  The  null  is  near  A  =  85°,  because  A  is  large  enough  that 
the  exponential  factor  must  be  considered.  At  greater  heights,  shown 
in  figures  3.  5,  3.  6,  etc.  ,  the  lobes  are  closer  together  as  expected. 

Figures  3.  7  and  3.  8  have  an  expanded  angle -of-elevation  scale 
to  show  the  detail  at  small  angles.  At  these  heights,  there  is  more  low- 
angle  frequency  dependence,  and  even  horizontal  polarise ati on  shows 
appreciable  amplitude  for  A  <  0.  For  0  <  A  <  5°,  the  propagation  factor  is 
greater  for  horizontal  polarization  than  for  vertical  polarization;  vertical 
polarization  is  better  at  negative  angles. 

The  propagation  factors  for  poor  land  are  shown  in  figure  4  for 
the  same  set  of  heights.  There  is  a  little  more  frequency  dependence 
at  low  heights,  and  the  difference  between  the  polarizations  is  not  so 
great.  At  medium  heights  the  behavior  of  the  propagation  factor  for 
vertical  polarization  is  complicated  by  the  psuedo-Brewster  angle  of 
the  ground  reflection  coefficient  [Jordan,  1950],  This  occurs  near 
A  =  6°  at  1  MHz  and  near  A  =  17°  at  1  5  and  2  5  MHz. 

For  greater  heights,  the  lobe  structure  is  similar  to  that  for  sea 
water  but  the  range  between  a  maximum  and  an  adjacent  minimum  is  not 
so  great.  Horizontal  polarization  is  slightly  better  thm  vertical  polari¬ 
zation  at  low  angles, 
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2.  7  The  Effect  of  a  Large  Irregularity 
So  far,  it  has  been  assumed  that  the  earth  is  smooth  and  homoge¬ 
neous.  All  parts  of  the  real  world  do  not  satisfy  this  requirement  so  the 
natural  question  is:  How  much  will  irregularities  affect  the  propagation 
factor?  A  complete  answer  to  this  question  is  not  available,  although  a 
number  of  special  cases  have  been  treated  in  the  literature.  In  general, 
it  is  assumed  that  surface  irregularities  very  small  compared  to  a  wave¬ 
length  can  be  accounted  for  by  using  ’'effective”  ground  constants. 

Structures  and  irregularities  comparable  to  a  wavelength  may  act  as 
scatterers,  and  large  irregularities  result  in  shadowing  and  diffraction. 

A  limited  attempt  to  compute  the  effect  of  a  certain  class  of  large 
irregularities  is  described  in  this  section. 

The  nature  of  the  irregularity,  a  large  box-like  ridge  or  mesa, 

is  illustrated  in  figure  5.  The  antenna  is  at  A,  d  kilometers  from  the 

H 

hill  of  height  h  .  The  solution  for  this  model,  given  by  Furutsu  [1965], 

is  a  triple  residue  series  with  terms  similar  to  those  in  (32).  This 

solution,  developed  for  VHF  tropospheric  paths,  is  not  well-adapted  to 

antenna  pattern  calculations  because  the  numerical  procedures  wilL  not 

work  except  for  very  small  elevation  angles.  However,  it  can  be  used 

to  determine  how  small  and  how  remote  an  irregularity  must  be  before 

the  propagation  factor  calculations  are  correct  at  low  angles. 

The  propagation  factor  was  calculated  for  lov  antennas  and  for 

elevation  angles  less  than  one  degree  for  hill  sizes  ranging  from  h  =  X/3 

H 

to  h^  =  10X  and  for  distances  d^  ranging  from  2.  5  km  to  10  km.  The  results 
were  divided  by  the  propagation  factor  for  a  spherical  earth  of  the  same 
conductivity  to  minimize  the  influence  of  ground  conductivity.  (This 
also  removed  the  polarization  dependence.  )  The  results  indicated  th^ 
following  empirical  rule:  For  the  small  angles  where  the  two-wave  theory 
cannot  be  used,  an  irregularity  of  the  kind  pictured  in  figure  5  will  cause 
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the  propagation  factors  to  be  in  error  by  1  dB  or  more  if  the  ratio  h  /d 

H  H 

exceeds  0.  01,  i.  e.  ,  if  the  elevation  angle  to  the  horizon  is  more  than 
one-half  degree. 

For  larger  angles,  no  calculations  could  be  made,  but  since  the 
two-wave  theory  is  valid  for  large  angles,  the  notion  of  requiring  a 
smooth  unobstructed  first  Fresnel  zone,  or  portion  thereof,  as  outlined 
by  Utlaut  [  1 962 ]  can  probably  be  used.  In  this  regard,  Anderson  [1964] 
has  found  that  the  Fresnel  zone  is  considerably  smaller  over  poorly 
conducting  ground  than  it  is  over  perfectly  conducting  ground. 

Much  more  research  on  the  effect  of  terrain  irregularities  is 

needed. 

3.  Antenna  Patterns 

In  the  preceding  section,  formulas  were  developed  for  propa¬ 
gation  factors  for  very  short  dipoles.  Most  HF  antennas  are  comparable 
to  a  wavelength,  and  it  is  necessary  to  integrate  along  their  length  to 
find  the  radiation  pattern.  This  was  done  for  several  HF  antennas 
located  above  a  flat  earth  by  Laitinen  [1949].  In  this  section,  the 
equations  for  vertical  and  horizontal  linear  antennas  over  curved  earth 
will  be  derived,  with  sample  patterns  showTi  for  X/  2  vertical  and  hori¬ 
zontal  dipoles  and  X  / 4  vertical  monopoles. 

3.  1  Vertical  Antenna  Pattern  Formulas 

A  thin  current-carrying  wire  of  length  L  =  h2  -  hi  is  located  on 
a  radial  of  the  earth  (h2  is  the  height  to  the  top  of  the  antenna;  hi  is  the 
height  to  the  bottom  of  the  antenna).  The  current  in  the  wire  is  assumed 
to  be  distributed  sinusoidally,  i.  e.  , 
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(42) 


1(h)  =  I  sink  (he  -  h)  , 
m 

where  I  is  the  current  maximum  and  h  =  r  -  a.  From  section  2,  the 

m  3 

radial  field  due  to  an  element  dh  is 


dE  —  E0  cos2A  I  sin[k(hg  ~h)]  F  (A,h)  dh, 
r  m  r 


(43) 


where  E0  =  -30  ik 


-ikD 

* 

D 


The  phase  path  to  T  is  a  function  of  position  along  the  antenna,  and 
this  is  accounted  for  by  multiplying  (43)  by  e  i.  e.  ,  the  phase 

is  referred  to  the  surface  of  the  earth.  Then 

h2 

ikhsinA  . 


E  —  E0  cos  A  I 


m 


i 


sin[k(h2  -h)]  F  (A,h)  dh.  (44) 


When  the  elevation  angle  A  is  large  enough  to  use  (5),  the  integral  can  be 
evaluated  in  closed  form,  and 


ikh2  sin  A 

E  —  E0  - - - -  I  f"  l+e 

r  cos  A  m  L 


-ik2h2  sinA 


R, 


_e-ikLsinAj.sinA  ginkL  q  _  e.ik2hl  sinAR^ 
+COS  kL.  (l  +  e-2ikhl  SinARv)}] 


ikh2  sinA  . 

-  E0  --0'8  A  Im[Fr(A,h2)-e'lkLsin  {i  sinAsinkLF0(A,hl) 

+  cos  kL  F^  ’  A,hi )  J-  I  .  (45) 

3.  The  formulas  in  this  section  are  not  complete  if  A  is  too  near  90°. 
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These  equations  show  that  if  the  elevation  angle  is  large  enough,  the 
field  is  proportional  to  the  propagation  factor  evaluated  at  the  top  of 
the  antenna  plus  the  propagation  factor  at  the  bottom  (with  the  phase 
shifted  to  account  for  the  difference  in  height).  This  reveals  the  basic 
nature  of  the  propagation  factors.  In  fact,  for  small  A, 


Fr(A.^)+e-iII-inArr(A.hl)-Fr(A,ii^). 


i.  e.  ,  the  antenna  pattern  of  a  vertical  half-wave  dipole  is  approximately 
proportional  to  the  propagation  factor  calculated  at  the  midpoint.  The 
accuracy  of  this  approximation  is  studied  later;  it  is  very  useful  for 
understanding  the  qualitative  behavior  of  the  antenna  pattern. 

For  small  A,  the  propagation  factor  is  itself  an  integral  and 
(44)  cannot  be  integrated  in  closed  form.  However,  the  integration  is 
easily  done  numerically  on  a  computer.  Patterns  of  the  two  types  of 
vertical  antennas  will  be  discussed  after  the  normalization  is  defined. 
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3.  2  Normalization  Relative  to  an  Isotropic  Source 

One  common  reference  antenna  is  the  (non-physically  realizable) 

isotropic  source  --an  antenna  that  radiates  equally  in  all  directions. 

The  E  field  of  an  isotropic  source  is  [Jordan,  1950] 
r 


-ikD 

(isotropic)  =  */3 OP  — g cos  A 

where  P  is  the  radiated  power  in  watts.  Let  R  ,  be  the  radiation 

rad 

resistance  of  the  antenna  being  compared.  The  current  I  is  found  from 

p  =  i2r  . 

rad 


and  the  pattern  of  a  vertical  antenna  relative  to  the 'isotropic  radiator  is 


E 

_ r _ 

E  (isotropic) 
r 


—  -  ik  cos 


dh. 


(49) 


The  impedance  of  a  half-wave  dipole  over  ground  is  the  sum  of 

its  free  space  self -impedance  plus  the  mutual  impedance  Zxg  between 

the  antenna  and  its  (imperfect)  image  in  the  ground.  According  to 

Jasik  [1961],  Z21  is  the  product  of  the  mutual  impedance  between  the 

antenna  and  its  perfect  image  ,  and  the  ground  reflection  coefficient  for 

vertical  incidence.  The  real  part  of  the  impedance  then  is  R  ,.  The 

rad 

radiation  resistance  for  a  vertical  half-wave  dipole  as  a  function  of 
height  of  the  midpoint  above  the  ground  is  shown  in  figure  6  for  the  types 
of  ground  in  this  report.  The  resistance  does  not  differ  much  from  the 
free  space  value  above  h~  I4A,  and  the  variation  is  greater  for  better 
conducting  ground. 
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3.3  Radiation  Patterns  of  Vertical  Half-wave  Dipoles 
Vertical  radiation  patterns  of  half-wave  dipoles  over  the  five 
types  of  ground  were  calculated  for  feed-point  heights  ranging  from 
.375  X  to  20  X  for  f  =  10  MHz  and  are  shown  in  figures  10-14.  The 
power  gain,  relative  to  an  isotropic  radiator,  in  the  direction  A  is 


G(A)  =  10  log10 


E  (A) 

r _ 

E^(isotropic) 


dB.  (50) 


This  is  the  function  plotted  in  the  figures. 

The  cha.racteristics  of  the  antenna  patterns  are  largely  deter¬ 
mined  by  the  behavior  of  the  ground  reflection  coefficient  Ry.  The 
amplitude  and  phase  of  Ry  for  the  frequency  and  ground  constants  of 
interest  are  shown  in  figure  8.  For  all  types  of  land  |Ry|  =  1  for  A  =  0°, 
and  decreases  to  a  minimum  as  A  increases  to  the  pseudo-Brewster 
angle  [Jordan,  1950].  (For  brevity,  the  pseudo-Brewster  angle  will 
hereafter  be  called  the  Brewster  angle.  )  The  Brewster  angle  is  less 
than  1°  for  sea  water  (at  10  MHz)  and  increases  with  decreasing  con¬ 
ductivity.  Above  the  Brewster  angle  \r\  recovers  and  approaches  its 
value  for  vertical  incidence,  which  is  greater  for  higher  conductivity. 

The  phase  of  the  reflection  coefficient  is  -180°  for  A  =  0  (R^  =  -1  when 
A  =  0),and,  except  for  the  curve  marked  ’’polar  ice”,  the  phase  increases 
slowly  until  it  reaches  the  vicinity  of  the  Brewster  angle.  There  it 
increases  quickly,  passing  through  -90°  at  the  Brewster  angle,  and 
then  flattens  out  to  its  value  for  vertical  incidence.  For  polar  ice,  the 
phase  decreases  at  first  to  -245°,  then  increases  quickly  to  about  -90°. 
This  different  behavior  of  the  phase  of  R  for  polar  ice,  together  with 
the  negligible  amplitude  of  R  above  A  =  30°,  explains  some  apparently 
anomalous  antenna  patterns  to  be  seen  later. 
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Now,  using  (48)  and  (1),  the  qualitative  behavior  of  the  antenna 
patterns  can  be  deduced.  First,  for  very  small  heights  h, 

Fr(h,  A)^l+Ry.  (51) 

When  A  is  less  than  the  Brewster  angle,  R^.  partially  cancels  the  1  (the 

reflected  wave  partially  cancels  the  direct  wave)  since  the  phase  of 

is  between  -180°  and  -90°.  The  extent  of  the  cancellation  depends  on 

|R  j  ,  so  in  this  region  |F  |  is  smaller  for  higher  conductivities. 

V  r 

Above  the  Brewster  angle,  the  reflected  wave  reinforces  the  direct 
wave,  and  the  field  exceeds  the  free  space  value. 

For  larger  heights,  the  phase  factor, 

-icp  -ik2h  sin  A  , 

is  important  in  determining  whether  the  reflected  ray  reinforces  or 
cancels  the  direct  ray.  Even  if  R^  =  -1,  the  reflected  ray  adds  if 
2kh  sin  A  =  (2n  -  1)  tt  for  n  =  1,  2  ,  ....  For  great  heights,  the  first 
lobe  of  the  antenna  pattern  is  near 

sin  A  —  •—  (52) 

if  A  is  below  the  Brewster  angle.  If  A  is  above  the  Brewster  angle,  the 
solution  of  (52)  is  near  a  null. 

This  general  qualitative  behavior  of  Fy  can  now  be  used  in  the 
interpretation  of  the  vertical  half-wave  dipole  patterns  in  figures  10-14. 
Figure  10  shows  the  vertical  radiation  patterns  in  dB  above  an  isotropic 
radiator  for  the  antenna  over  sea  water  for  various  heights  of  the  feed- 
point,  ^(hi+hs).  Figure  10.1  is  for  the  lowest  height,  §-\;  figures  10.2 
10.3,  10.4,  10.5,  10.6,  1 0.  7,  and  1 0.  8  are  for  feed  point  heights  of  j,  J-,  1, 
2,  5,  10,  and20\,  respectively.  The  frequency  is  10  MHz.  The 
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calculations  were  made  with  both  the  flat  and  spherical  earth  formulas 
in  the  region  0  <  A  <  5°,  and  the  flat  earth  values  are  indicated  by  a 
diamond.  In  figure  10.  1,  the  maximum  of  the  pattern  is  about  7  dB  at 
A  =  6°.  The  Brewster  angle  for  sea  water  is  at  A  <  1 0 ,  so  the  reflected 
wave  is  reinforcing  the  direct  wave  and  the  gain  is  well  above  the  2.  15 
dB  gain  of  a  half-wave  dipole  in  free  space.  It  does  not  reach  the  8.  17 
dB  gain  expected  over  perfect  ground  because  approaches  unity  only 
for  large  A. 

Note  that  there  is  appreciable  power  at  A  =  0  --  about  equal  to 
an  isotropic  --  and  even  some  for  negative  elevation  angles. 

The  picture  is  very  different  for  the  same  height  antenna 
over  good  ground,  shown  in  figure  11.  1.  Here  the  Brewster  angle  is 
near  A  =  14°,  so  at  lower  angles  the  reflected  wave  is  partially  cancelling 
the  direct  wave.  The  pattern  reaches  its  maximum  of  0  dB  near  14c  where 
|RVI  is  smallest  and  there  is  least  cancellation.  One  might  expect  the 
gain  to  be  greater  than  the  free  space  value  for  A  >  14°,  since  the  phase 
of  Ry  is  such  that  j  1  +  R^|  >  1.  However,  A  ismow  large  enough  that 
the  phase  path  difference  is  now  significant.  For  example,  at  A  =  20°, 
cp  =  2kh  sin  A  ^  1. 61  radians  —  92°.  The  propagation  factor  is  then  about 
1  +  |Ryl  e  1  Jits  amplitude  is  less  than  one,  and  the  gain  is  there¬ 
fore  less  than  in  free  space.  This  illustrates  a  simple  rule:  When  both 
A  and  h  are  very  small,  the  behavior  of  1  +  R^  qualitatively  describes 
the  antenna  pattern  fairly  well,  but  if  either  A  or  h  is  large  (more 

precisely,  if  h  sin  A  is  greater  than,  say,  •—  )  the  phase  path  difference 

-icp 

factor,  e  ,  must  also  be  considered. 

The  antenna  at  the  same  height  over  poor  ground  has  the  same 
general  behavior  as  shown  in  figure  12.  1.  The  maximum  gain  is 
slightly  greater  than  for  good  ground  since  |Ry|  is  smaller  and  less 


27 


of  the  direct  wave  is  cancelled.  The  maximum  gain  occurs  at  a  higher 
angle,  A  —  19°,  since  the  Brewster  angle  is  at  a  larger  A. 

This  trend  is  continued  as  shown  in  figure  13.1  where  the 
conductivity  is  very  low.  The  maximum  gain,  about  3  dB,  now  exceeds 
the  gain  in  free  space  because  of  the  different  behavior  of  the  phase  of 
Ry  that  was  mentioned  earlier.  The  pattern  is  very  near  to  that  for 
free  space  above  A  “  30°  because  |Ry|  is  so  small. 

The  Brewster  angle  for  fresh  water  is  near  7°,  so  the  gain  is 
less  than  it  is  in  free  space  for  A  <  7*  in  figure  14.  1.  The  maximum  is 
about  3.  5  dB  above  an  isotropic  radiator  at  A  =  12°. 

For  a  feedpoint  height  of  X/2,  (52)  predicts  a  null  at  A  =  30°, 
which  can  be  seen  in  figure  10,2  There  is  a  minimum  here  in  figures 
11.  2  and  14.  2,  although  it  is  not  so  deep  because  |Ry|  is  smaller  and 
less  of  the  direct  wave  is  cancelled.  In  figures  12.  2  and  13.2  the  effect 
can  hardly  be  seen  because  |Ry|  is  so  small. 

The  pattern  of  a  vertical  dipole  fed  f-  X  above  the  ground  has  the 
same  general  characteristics  as  those  for  the  two  lower  heights.  In  the 
conclusions,  these  three  cases  are  grouped  as  "low  heights". 

Similarly,  heights  of  X  and  2X  are  grouped  as  "medium  heights" 
and  only  the  one  wavelength  case  will  be  discussed.  If  A  is  above  the 
Brewster  angle,  (52)  predicts  a  null  near  14^°  which  can  be  seen  in 
figures  10.  4  and  14.  4  for  sea  water  and  fresh  water,  respectively.  The 
pattern  for  sea  water  has  a  good  lobe  at  3°  which  is  above  the  Brewster 
angle.  The  fresh  water  pattern  has  a  relative  maximum  near  7R,  where 
the  gain  is  about  that  in  free  space  because  this  is  also  the  location  of 
the  Brewster  angle.  The  Brewster  angle  for  good  ground  is  about  12°, 
so  the  pattern  for  this  case  (figure  11.4)  does  not  have  a  relative 
minimum  there  and  is  almost  equal  to  the  free  space  value.  As  expected, 
the  pattern  for  poor  ground  is  between  those  for  good  ground  and  polar  ice. 


28 


The  patterns  for  feed  point  heights  of  5,  10,  and  2  OX  have  similar 

characteristics  and  are  grouped  as  "high  heights".  The  highest,  h  =  20X 

shown  in  figures  10.  8,  11.8,  etc.  ,  will  be  discussed  as  typical.  The 

-icp 

phase  of  the  reflected  ray,  e  R^,  is  now  dominated  by  the  phase  path 
difference  cp.  Thus,  all  cases  except  sea  water  have  a  lobe  near  0.7° 
predicted  by  (52),  with  a  gain  of  about  7  dB.  The  Brewster  angle  for 
sea  water  is  near  0.  6®  and  this  moves  the  position  of  the  first  lobe 
above  1°.  At  these  high  heights,  the  pattern  oscillates  around  the  free 
space  value,  and  the  envelope  of  the  pattern  is  also  of  interest.  The 
Brewster  angle  is  indicated  by  a  narrowing  of  the  envelope,  since  the 
interference  is  proportional  to  |Ry|.  Figures  12.  8  and  14.8  are 
especially  good  examples. 

The  success  in  explaining  the  behavior  of  the  antenna  patterns 
using  the  propagation  factor  for  the  midpoint  of  the  vertical  dipole  suggests 
that  this  may  be  a  very  good  approximation.  The  approximation  (48)  was 
tested  numerically  for  small  take-off  angles,  with  the  additional  assumption 
that  the  radiation  resistance  of  the  antenna  is  approximately  its  free 
space  value.  The  formula  is 


G(A)  S=  2.  15  +  20  logic  |F  <A,-^4^))|.  (53) 

r  ^ 

/ 

For  0  <  A  <  5°,  the  error  was  less  than  1  dB,  and  less  than  2  dB  for 
negative  elevation  angles.  The  approximation  was  not  tested  above 
A  =  5°  since  low  angles  are  of  most  interest  here,  but  the  approximation 
is  probably  adequate  up  to  A  =•  30°. 
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3.4  Radiation  Patterns  of  Quarter -Wave  Monopoles 


Equation  (47)  gives  the  field  of  a  quarter -wave  grounded  antenna 
over  flat  earth,  and  (44)  can  be  integrated  numerically  over  curved 
earth.  We  assume  a  perfect  counterpoise  or  ground  screen  large 
enough  to  sustain  the  sinusoidal  current  distribution,  but  too  small  to 
affect  the  pattern.  The  radiation  resistance  is  36.6  ohms,  so  that  the 
pattern  relative  to  an  isotropic  radiator  is 

E  ft/4)  h2 

- -  -0.  905  ik  cos  A  |  . . .  dh 

r(isotropic)  hx 

-0.  905  i  i  k  h2  sin  A 

~ - : — e 

cos  A 

r  .  TT  .  . 

x|Fr(h2)  -isinAe'V111  Fgtfi!)}  .  (54) 


The  second  line  is  for  flat  earth.  The  term  involving  Fg  can  be  inter¬ 
preted  as  the  field  reradiated  by  the  horizontal  ground  system  since  it 

is  the  same  as  the  propagation  factor  for  the  vertical  field  off  the  end 
of  a  horizontal  electric  dipole  [Laitenen,  1949]. 

Patterns  were  calculated  for  a  frequency  of  10  MHz  for  all  ground 
types  (including  perfectly  conducting  ground)  for  two  heights  of  the  feed- 
point:  hi  =0,  figure  15,  and  hi  =  X/4,  figure  16.  For  these  low  heights 
the  pattern  is  roughly  proportional  to  1  +  for  small  A.  Thus,  below 
the  Brewster  angle,  the  reflected  ray  tends  to  cancel  the  direct  ray; 
above  the  Brewster  angle  it  reinforces  the  direct  ray.  The  gain  for  sea 
water  is  almost  as  good  as  for  perfect  ground  but  comes  at  a  higher 
angle.  Then,  for  successively  higher  Brewster  angles,  come  fresh 
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water,  good  ground,  and  poor  ground.  As  for  the  half-wave  dipole  at 

low  heights,  the  gain  for  good  and  poor  ground  never  reaches  the  free 

space  value  because,  by  the  time  the  Brewster  angle  is  exceeded,  A  is 

-icp 

large  enough  that  the  factor  e  causes  partial  cancellation.  The  pattern 
for  polar  ice  is  very  similar  to  that  for  poor  ground. 

The  patterns  for  a  quarter -wave  monopole  with  its  counterpoise 
a  quarter -wave  above  the  ground  are  shown  in  figure  16.  These  have  a 
relative  minimum  between  45°  and  50°.  Otherwise,  they  have  the  same 
features  seen  in  figure  15.  Over  sea  water,  the  power  is  only  3  dB 
below  an  isotropic  radiator  at  A  =  0°,  with  some  power  at  negative  take¬ 
off  angles. 

3.  5  Vertical  Patterns  Broadside  to  a  Horizontal 
Half-Wave  Dipole 

When  the  linear  antenna  is  horizontal,  there  is  no  variation  of  the 
current  with  h,  and  the  vertical  radiation  pattern  is  given  exactly  by  the 
free  space  pattern  multiplied  by  the  propagation  factor  for  horizontal 
polarization  [Wait*  1956],  The  field  broadside  to  a  horizontal  dipole 
(in  the  plane  perpendicular  to  the  dipole  through  its  midpoint)  relative 
to  an  isotropic  radiator  is 


y.mt  . (h 

E  (isotropic)  rz - cp  A 


cp 


VR 


rad 


(55) 


The  power  gain  varies  inversely  with  the  radiation  resistance 
which  is  shown  in  figure  7  for  a  horizontal  half -wave  dipole  over 
various  types  of  ground.  The  variation  from  the  free  space  value  is 
greater  than  that  for  the  vertical  antenna. 


s 
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To  help  understand  the  patterns,  the  behavior  of  the  reflection 

coefficient  for  horizontal  polarization,  R  ,  is  reviewed.  The  amplitude 

h 

and  phase  of  R,  are  shown  in  figure  9  as  a  function  of  A  for  a  frequency 
h 

of  10  MHz  and  for  the  ground  types  of  table  1.  There  is  no  psuedo-Brewster 

angle  for  this  polarization,  and  the  most  important  feature  is  that  R  is  *  ' 

h 

not  too  different  from  -1  for  any  take-off  angle  less  than  2.0°,  except 
for  polar  ice.  Therefore,  if  the  phase  path  difference  cp  is  small, 
the  reflected  ray  nearly  cancels  the  direct  ray.  The  extent  of  the 
cancellation  will  depend  on  Ir^  so  that  the  antenna  will  ha^e  a  better 
low  angle  pattern  over  lower  conductivity.  The  obvious  way  to  increase 
the  phase  path  difference  and  thus  improve  the  low  angle  pattern  is  to 
increase  h^  (raise  the  dipole). 

The  calculated  vertical  radiation  patterns  of  a  horizontal  half¬ 
wave  dipole  located  0.  05 A.  above  the  five  ground  types  are  shown  in 
figure  17.  1.  As  expected  from  the  above  discussion  of  reflection 
coefficients,  very  little  power  is  radiated  at  low  angles.  For  all  ground 
types  except  polar  ice,  the  power  is  more  than  10  dB  below  an  isotropic 
radiator  if  A  <  10°.  The  radiation  is  less  over  the  highly  conducting 

surfaces  (such  as  water)  because  R  is  nearer  to  -1  and  the  reflected 

h 

ray  more  nearly  cancels  the  direct  ray. 

As  the  antenna  is  raised  to  0.  IX,  figure  17.2,  and  the~  to  0.  25X, 
figure  17.3,  the  radiation  increases  but  has  the  same  general  character¬ 
istics.  The  better  conducting  models  now  have  more  radiation  at  high 
angles  because  the  phase  path  difference  is  large  enough  that  the 
reflected  ray  reinforces  the  direct  ray.  This  trend  continues  at 
h^  =  0.  375X  and  X/2  as  shown  in  figures  17.  4  and  17.  5.  N^w,  cp  is 
large  enough  (for  some  values  of  A)  that  the  antenna  gain  over  an  isotropic 
is  more  than  8  dB. 
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As  the  antenna  is  raised  higher  and  higher,  as  shown  in  figures  17.  6 

17.  10,  a  standing  wave  pattern  is  seen  with  lobes  at  elevation  angles  A  , 

»  n 

such  that 

sin  A  —  (2n  -1)  —  . 
n  4h 

The  gain  is  best  for  sea  water.  The  pattern  for  an  antenna  over  polar 
ice  is  a  little  different  from  the  others  for  large  A  because  |RjJ  is  very- 
small  for  this  case. 

3.6  Comparison  of  Low-Angle  Patterns  of  the 
Three  Types  of  Antennas 

All  the  antenna  patterns  show  power  gain  relative  to  an  isotropic 
radiator  so  they  can  be  compared  directly.  In  this  section,  their  relative 
effectiveness  at  low  elevation  angles  is  discussed. 

The  half-wave  vertical  dipole  fed  §-\  above  the  ground  radiates 
better  at  low  angles  than  the  quarter -wave  monopole.  In  a  sense,  this 
comparison  is  not  "fair",  since  the  discussion  above  showed  that  height 
is  a  great  advantage  in  attaining  a  good  low-angle  pattern,  and  this  half¬ 
wave  dipole  is  higher  than  the  quarter -wave  monopole.  The  quarter -wave 
monopole  has  a  better  low-angle  pattern  than  the  horizontal  dipole  if  the 
dipole  is  less  than  a  half-wave  above  the  surface,  except  for  polar  ice. 
Over  sea  water,  the  monopole  is  better  than  the  horizontal  dipole  at  very 
low'  angles  (<  5°  )  unless  the  dipole  is  more  than  5^  high,  although  the 
maximum  gain  of  the  dipole  is  greater. 

Now'  compare  the  low-angle  patterns  of  horizontal  and  vertical 
dipoles  fed  at  the  same  height.  At  the  lowest  height  for  which  the 
vertical  dipole  was  computed,  h  =  f-\,  the  vertical  dipole  has  a  better 
low-angle  pattern  than  the  horizontal  over  all  types  of  ground  except 
polar  ice.  At  the  other  low  and  medium  heights,  the  vertical  dipole  is 
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better  over  water,  and  the  horizontal  is  better  over  land.  For  high 
antennas,  the  horizontal  dipole  has  a  slightly  better  pattern  at  very  low 

•x 

angles. 

3.  7  Remark  on  Frequency  Dependence 
The  antenna  patterns  in  this  report  were  calculated  for  a  frequency 
of  10  MHz,  so  it  may  seem  that  the  results  are  of  limited  application. 
However,  the  graphs  of  the  propagation  factor  in  section  2  show  that 
there  is  only  a  weak  dependence  on  frequency,  and  these  same  graphs 
can  be  used  to  estimate  the  frequency  dependence  of  the  antenna  patterns. 
When  doing  this,  remember  that  all  heights  have  been  measured  in  wave¬ 
lengths.  An  antenna  at  some  fixed  physical  height  would  have  more 
frequency  dependence  since  its  height  in  wavelengths  changes  with 
frequency. 

4.  Conclusions 

Propagation  factors  are  basic  to  the  calculation  of  vertical 

radiation  patterns  of  high  frequency  antennas  over  finitely  conducting 

earth.  For  elevation  angles  greater  than  one  or  two  degrees,  they  are 
-iCp 

given  by  1  +  e  R  where  R  is  the  appropriate  ground  reflection  coefficient 
and  cp  is  the  difference  between  the  phase  paths  of  the  direct  and  reflected 
rays.  For  smaller  elevation  angles  (including  negative  elevation  angles), 
the  full-wave  solution  for  diffraction  by  a  spherical  earth  must  be  used. 

The  calculations  in  this  report  show  the  effects  of  frequency, 
conductivity,  polarization,  and  height  on  the  propagation  factor.  Over 
sea  water  there  is  very  little  frequency  dependence  for  positive  elevation 
angles.  There  is  more  variation  with  frequency  over  land,  especially 
for  vertical  polarization.  The  conductivity  has  a  much  greater  effect 
on  vertical  polarization  than  on  horizontal  polarization.  In  general,  for 
low  heights,  the  amplitude  of  the  vertically  polarized  propagation  factor 
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grows  as  the  conductivity  increases  and  the  amplitude  of  the  horizontally 
polarized  propagation  factor  decreases  with  increasing  conductivity.  At 
great  heights,  the  propagation  factors  for  both  polarizations  improve 
with  increasing  conductivity. 

The  vertical  radiation  pattern  broadside  to  a  horizontal  dipole  is 
proportional  to  its  propagation  factor.  This  is  not  exactly  true  for  a 
vertical  dipole,  but  the  error  is  less  than  1  dB  for  positive  elevation 
angles  if  the  propagation  factor  for  the  mid-point  of  the  antenna  is  used. 
Patterns  were  calculated  for  vertical  and  horizontal  half-wave  dipoles 
and  vertical  quarter-wave  monopoles.  For  low  heights,  the  vertical 
dipole  has  the  best  low-angle  radiation  pattern;  for  medium  heights, 
the  vertical  dipole  has  the  best  low-angle  pattern  over  water,  and  the 
horizontal  dipole  is  better  over  land.  The  horizontal  dipole  is  better 
for  all  types  of  ground  at  great  heights.  The  best  low-angle  gain, 
relative  to  an  isotropic  radiator  is:  for  low  antennas,  about  7  dB  over 
sea  water  and  about  0  dB  over  land;  for  medium  heights,  about  6  dB; 
and  for  high  antennas,  about  8  dB. 

Further  research  is  required  to  predict  the  quantitative  effect 
of  irregularities  on  low-angle  radiation  patterns. 
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opagation  factors  over  sea  water, 
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RADIATION  RESISTANCE 
OF  A  VERTICAL  X/2  DIPOLE  OVER  GROUND 


6.  Radiation  resistance  of  a  vertical  half-wave  dipole  over  various  types  of 
ground  as  a  function  of  height  of  the  midpoint. 
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RADIATION  RESISTANCE 
OF  A  HORIZONTAL  X/2  DIPOLE  OVER  GROUND 
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POLAR  ICE 


take-off  angle,  degrees 

Phase  of  the  ground  reflection  coefficient  for  vertical  polarization, 
R  ,  for  various  types  of  ground.  The  frequency  is  10  MHz. 


SEA  WATER 


PHASF 


POLAR  ICE 


TAKE-OFF  ANGLE,  DEGREES 

9.  2.  Phase  of  the  ground  reflection  coefficient  for  horizontal  polarization 

R,  ,  for  various  types  of  ground.  The  frequency  is  10  MHz. 
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Vertical  Halfwave  Dipole 
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10,  1.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located 
above  sea  water.  The  frequency  is  10  MHz. 
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IEI  <RELATIVE  TO  ISOTROPIC)  IN  OB 

10.4.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  X 
above  sea  water.  The  frequency  is  10  MHz. 
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10.  5.  Vertical  radiation  pattern  of  a  vertical  half-wave  ai; 
above  sea  water.  The  frequency  is  10  MHz. 
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10.6.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  5X 


above  sea  water.  The  frequency  is  10  MHz. 
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10.7.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  10X 
above  sea  water.  The  frequency  is  10  MHz. 
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90-  Frequency  =  10  MHz 

Vertical  Halfwave  Dipole 
Feed  Point  Height  =  600.  00  Meters 
Sea  Water 
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I  E  I  RELATIVE  TO  ISOTRJPIO  IN  DB 

10.  8.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  20 X 
above  sea  water.  The  frequency  is  10  MHz. 
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11.1.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  f*  X 
above  good  ground.  The  frequency  is  10  MHz. 
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11.2.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  k/Z 
above  good  ground.  The  frequency  is  10  MHz. 
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11.5.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  Z\ 
above  good  ground.  The  frequency  is  10  MHz. 
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11.6,  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  5X 
above  good  ground.  The  frequency  is  10  MHz. 
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11.7.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  10X. 
above  good  ground.  The  frequency  is  10  MHz. 
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IEI  <RELATIVE  TO  ISOTROPIO  IN  DB 

Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  loca.jd  20 X 
above  good  ground.  The  frequency  is  10  MHz. 
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IEI  <RELATIVE  TO  ISOTROPIO  IN  DB 

12.  3.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  f-  X 
above  poor  ground.  The  frequency  is  10  MHz. 
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12.  5.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  2X 
above  poor  ground.  The  frequency  is  10  MHz. 
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12.6.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  5X 
above  pc  or  ground.  The  frequency  is  10  MHz. 
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Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  f-A. 
above  polar  ice.  The  frequency  is  10  MHz. 
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13.  2.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  X/2 
above  polar  ice.  The  frequency  is  10  MHz. 
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13.4. 


Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  X 
above  polar  ice.  The  frequency  is  10  MHz. 
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13.6.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  5X 
above  polar  ice.  The  frequency  is  10  MHz. 
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Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  10X 
above  polar  ice.  The  frequency  is  10  MHz. 
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14.  5.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  2  A. 
above  fresh  water.  The  frequency  is  10  MHz. 
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14.  7.  Vertical  radiation  pattern  of  a  vertical  half-wave  dipole  located  10\ 
above  fresh  water.  The  frequency  is  10  MHz. 
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I  E  t  RELATIVE  TO  ISOTROPIC)  I  il  DD 

Vertical  radiation  pattern  of  a  quarter -wave  grounded  antenna  with  perfect 
counterpoise  on  various  kinds  of  ground.  The  frequency  is  10  MHz. 
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I  E  I  <RELATIVE  TO  ISOT..OPIO  IN  DB 
Vertical  radiation  pattern  of  a  quarter -wave  antenna  with  perfect 
counterpoise  located  X/4  above  various  types  of  ground.  The  frequency 
is  10  MHz. 
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17.  1.  Vertical  radiation  pattern  of  a  horizontal  half-wave  dipole  0.  05X  above 
various  kinds  of  ground.  The  frequency  is  10  MHz. 
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Vertical  radiation  pattern  of  a  horizontal  half-wave  dipole  Z\  above 
various  kinds  of  ground.  The  frequency  is  10  MHz. 
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17.9.  Vertical  radiation  pattern  of  a  horizontal  half-wave  dipole  5 A.  above 

various  kinds  of  ground.  The  frequency  is  10  MHz. 


Frequency  =  10  MHz 
Horizontal  Dipole 
Feed  Point  Height  =  600.  00  Meters 
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Vertical  radiation  pattern  of  a  horizontal  half-wave  dipole  20 X  above 
various  kinds  of  ground.  The  frequency  is  10  MHz. 
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